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Abstract 

The effective Hamiltonian, as obtained from applying the Hamiltonian flow equa- 
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1 Introduction 

Light-front frame is believed to be a useful tool for solving bound state problems in QCD. 
Generally, bound state calculations in the field theory have two sources of complexity - 
they are relativistic and of many-body type. The method of flow equations copes with 
both of them, at least to the definite order in perturbation theory. One transforms the 
Hamiltonian to eliminate interactions changing particle number, reducing thus the bound 
state problem to a few-body problem. Simultaneously utraviolet divergencies occur, origi- 
nating from the high-energy region. To complete renormalization one uses either coupling 
coherence or fixes counterterms to provide finite values for physical observables and to 
retain symmetries violated by the procedure 0. Both type of fiow equations of renormal- 
ization type and for the new (particle number conserving) interactions appear together. 
This program can be fulfilled in perturbation theory expansion. As a result, the bound 
state problem is approximated by a set of renormalized, effective interactions that do not 
change particle number. 

The sensitive tool in the light-front frame to check how accurately one desribes bound 
states by these effective interactions is to measure the violation of rotational symmetry. 
This symmetry is linked to a dynamical operator on the light-front, since rotations are 
dynamical, i.e. depend on the interaction. The symmetry may be spoiled in two steps: 
first, regularization and renormalization; and second, reduction to the effective few-body 
interactions with particle number conservation. The nonperturbative renormalization flow 
is of crucial importance for QCD p|, ||, |[ , but one can disregard this point in QED bound 
state calculations, that disentangles the two problems mentioned above. To the leading 
order the results for QED are obtained in where positronium system is described 
approximately by the effective electron-positron interaction. In the nonrelativistic limit 
the results for positronium spectrum agree with the results of covariant calculations. 

There are at least two other alternative approaches to solve for bound states in the 
light-front dynamics, the scheme of similarity renormalization of Glazek and Wilson 
and the method of iterated resolvents of Pauli [§. In both schemes calculations of the 
effective electron-positron interaction are performed and the question of rotational invari- 
ance for positronium spectrum is investigated. 

Calculations done so far in the similarity renormalization scheme use the nonrelativis- 
tic limit to flnd corresponding eigenvalues in the bound state perturbation theory [^. 
Analytical calculations are performed there for ground state: ground triplet levels are 
degenerate, indicating that rotational symmetry is restored [^. Performing similarity 
renormalization one eliminates high-energy modes and absorbs relativistic effects into an 
effective band-diagonal Hamiltonian, which describes bound state creation at nonrela- 
tivistic energy scales. It is a well working scheme for such systems as positronium 
therefore nonrelativistic approximations done in this approach to extract eigenvalues from 
effective Hamiltonian are quite natural there. In general, it is not always the case. In 
fact rotational symmetry becomes kinematic one, like light-front boost in the nonrela- 
tivistic limit, i.e. total momentum and its projection can be considered approximately as 
quantum numbers, that makes simpler to trace rotational invariance in these calculations. 

In the method of iterated resolvents an effective electron-positron potential is obtained 
and exact numerical solution of positronium bound state equation with the given poten- 
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tial is done [|^, 0. Degenerate multiplets for ground as well as for exited states are 
obtained |rD|. Q It is convenient to perform relativistic calculations in the light-front 
frame, which effectively has nonrelativistic kinematics. 

In the present work we perform relativistic few-body calculations for positronium 
spectrum numerically in the spirit of the work Trittmann et.al. |jlO[ (and using numerical 
code Ijllll), based on the effective electron-positron Hamiltonian obtained by the flow 
equations |jl5[. Effective interaction was derived there for different cutoff functions. The 
requirement of block-diagonalization of the Hamiltonian determines the generator only up 
to a unitary transformation of the blocks; this explains why the effective interaction may 
depend on the cutoff function. The question we investigate is to what extent rotational 
invariance is violated on the level of positronium spectrum and how does it depend on 
the choice of the cutoff function. We are not able to trace rotational invariance during 
the calculations, since it is dynamical operator. This is an excellent test for the method 
of flow equations itself and the control of approximations done during the calculations. 



2 Formulation of the problem 

We address to solve a light-front Hamiltonian bound state equation 

H\^P) = E\^) (1) 

for positronium. Using flow equations we transform the QED Hamiltonian H to a. block- 
diagonal effective Hamiltonian, which reduces positronium problem to a bound state 
problem in the electron-positron sector. The effective Hamiltonian for an electron and a 
positron is 

-f^efF = Hq + UeS (2) 

where Hq is the kinetic energy, and f/gfr includes effective interactions generated by the 
flow equations in the second order in coupling constant. The integral bound state equation 
is written 

-E'(Pi,P2; Ai,A2|^) = {Ep^+Ep^){pi,p2;Xi,X2\^) 

+ J2 Jd%d^P2{pi,P2; Ai, A2|f/effb'i,P2; A'l, X2){p'i,P2. a;, A'21^) 

(3) 

where the effective Hamiltonian pickes out from the positronium wave function the 
lowest ee-component {pi,P2] Ai, A2IV') with pi, Aj being the light-front three-momenta and 
helicities, respectively, carried by an electron {i = 1), and a positron {i = 2). The primed 
quantities refer to the initial state, the unprimed ones to the final state. The effective 
interaction {pi,p2', Ai, X2\Ues\Pi,P2] A'l, A2) = UesS{pi + P2 — Pi — P2) will be specified be- 
low. In order to deduce a Lorentz invariant energy we consider the bound state equation 
written for operator P^P^, corresponding to the invariant mass-squared on the light- 
front, rather than for the light-front Hamiltonian operator H = P^ . The light-front 

^ Numerical solution of positronium bound state problem in the light-front frame can be found also 

in 
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integral equation, 



Ai,A2| 



{x, K±;Xi, A2 



+ dx'd^K'^{x,Kr,>^l,>^2\Ves\x',K'^]X[,X2){x',K'^]X[,X2\ 



(4) 



is independent of the total momentum P"*" and P^. We introduced Ves = P'^'^Ucs- In that 
equation only intrinsic transversal momenta k± and longitudinal momentum fractions x = 
Pi/P~^ appear {p^ = {xP~^,xP± + i^±,Pi))- Its spectrum is thus manifestly independent 
of the kinematical state of the bound system, particularly of P"*" and Pj_, which reflects 
on the boost invariance peculiar to the light-front form |^. The integration domain D is 
restricted by the covariant cutoff condition of Brodsky and Lepage M, 



m 



xil 



X] 



(5) 



which allows for states having a kinetic energy below the bare cutoff A. The effective 
interaction between electron and positron, being a kernel in the integral equation (H), is 



generated by the flow equations |T5 

v;fr = 



a 



4vr2 

a 
A-n 



MM 



M2 




with the generator of unitary transformation 

1 fd\nf{D;iy 



D 



dl 



9(1)- 



(6) 



(7) 



where g{l) is the coupling constant as a function of flow parameter I, and f{D; I) is the 
cutoff function specifled below. In Eq. @ subscript ex refers to the exchange part, and 
an to the annihilation part. The null vector t]'^ has components [tj'^ ,'yfj_,V~) = (0,0,2) 
and is specific to the light-front calculations. The light-front metric tensor is denoted by 
Qfj^iy. The current-current tensors in the two channels are 

{u{pu X,)ru{p[, X[)) {v{p'2, X'2)rv{p2, A2)) 



\Jxx'{l — x)(l — x') 

iu{pu x,)rv {P2, A2)) {v{p'2, y2)ru{p[, a;)) 

xx'{l — — x') 



(8) 



where the fermion momenta were defined after Eq. (|^). The remaining definitions are as 
follows. The energy differences along the electron and the positron line. 



P2 



-Pi 

P>2 



(P'l - Pi 
{P2 



P2) 



(9) 
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respectively, have a simple relation to the (Feynman-) 4-momentum transfers along the 
two lines 

Ql = -ip[-pif = -q^D,, 

Ql = -{P2~p'2? = -q^D,- (10) 

Since the Feynman-momentum transfer Q is more physical quantity than the energy 
difference, we will use the former as far as possible. In fact, in our formulae we make use 
of the mean-square momentum transfer and the mean-square difference, 

SQ' = ^(Qe-a-)=-Ype-/^e-), (H) 

respectively. The dependence of the effective interaction Eq. (|^) on the cutoff function 
f{D] I) is carried by the factor 



oo 



0(D„ D,) = - I dl' ^i^£^f{D,; I') ^ Bee- , (12) 







which is asymmetric in the arguments but which satisfies 

Q{De, De) + Q{D,, De) = Qee + Qee = 1 • (13) 

The latter combination obeys 



Gee Oee ^ ( 

Ql Ql QlQl \ Q' ^ 

Gee- Q-ee SQ^ / 



Ql Ql QIQIK SQ 



(Gee- - Gee) • (M) 



that we use further. 

For the annihilation term we define the energy differences as 

Da = P'l' + P'2' - (P'l + P2)~ 

Db = Pi +P2 - {P1+P2Y ■ (15) 

They are related to the 4-momentum of the photon and to the free invariant mass- 
squares of the initial and final states 

Ml = {p[+p',r=p^D^, 

Ml = {pi+P2?=P^D,, (16) 
as well as to their mean and difference 

M' = \{Ml + Ml) = ^-^{Da + D,), 

SM' = \{Ml-Ml) = ^-{D^-D,), (17) 
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respectively. 

Effective interaction includes two different Lorentz structures: (7^1/ part insures 
the Bohr spectrum and is responsible for the spin splittings; rjfj^rji, term is diagonal in 
spin space and vanishes for real processes, i.e. on mass shell with SQ^ = 0, making the 
effective interaction to coincide with the Tamm-Dancoff approximation ||15|. The explicit 
x-dependence in the denominator of Eq. (||) looks like the only remnant of the light-front 
formulation; all other quantities are Lorentz scalars. One can absorb this dependence by 
redefining the wave function in the integral equation Eq. (^). We introduce instead of 
Jacobi momentum (x, k,±) the three momentum in the center of mass frame p = {pz, k,^) 
as follows 

2 ( ~'~ y/p'^ + m^) ' ^ 
where the Jacobian of this transformation dx / dp^ is 

2(p2+^2)3/2 ^ E ' ^ ' 

and it holds in this frame 



4 p ^ + ivP- 



x{l — x) 

E = ^Jp + w?. (20) 

The connection between 'old' and 'new' wave functions and the interaction matrix ele- 
ments are 



xil — x) 



(X, Kj_|Vefr|x , = , (21) 

\ix{\ — x)x'{\ — x') 



Integral equation 



W {p-KM^') = ^{p' + m^) {p;X,,X2m 

d^p' 



is written in a rotationally covariant form. 



E (P1 ^^iK'frlP a;, A^) {p'; a;, A'^I^') (22) 



3 Rotational invariance 

Integral equation (^) has rotationally covariant but still not rotationally invariant form 
because of the interaction kernel Veff, written in the light-front frame. Let us extract the 
part of interaction, which has manifestly rotational symmetry. 
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Quite generally factor Eq. (|l^) is a function of the ratio of its two arguments. 
Therefore, making use of Eq. (|T^), the effective interaction Eq. is given as 

Kff = -£^{YY)e.B,, - £^{YY)anC,. , (23) 

where the exchange part is defined 

m = Qe,-e-ee , e = ^, (24) 




and the annihilation part 



C„ 




M2 \ M2 p+2 ; 1 — /32 
X{P) = Qab-Qba , P=^- (25) 

The terms in the effective interaction proportional to "(9(0 5 xiP) depend explicitly on the 
choice of cut-off function and arise from /-ordering of the generator in the operator of 
unitary transform |T5|. Explicit form of the effective interaction with different cut-off 
functions is given in Appendix A. 

Define energy denominators in equation (^3]). Due to the three- momentum conserva- 
tion on the light-front, P1+P2 = Pi +P2 ioT longitudinal and transversal components, one 
has De — Dg = Da — -Db, where the energy denominators in both channels are given 
in Eq. (^ and Eq. (|15|). Therefore 

^ j _ (26) 

where - the total energy difference between initial and final states shows the 'off- 
shellness' of process. 

Using the parametrization Eq. (p!8[) one has for the energy denominators 

\Ma mJ 
Ml = 4(p'2 + m2) 

Ml = 4{p^ + m^), (27) 

where q = p'—p = [q^, q±) is the three- momentum transfer of the photon, and the relations 
between mean-squared and difference momenta and corresponding energy differences are 
given in Eq. ([Tl|) and Eq. ([T^ for exchange and annihilation channels, respectively. 

The second term in Eq. (p5|) is obviously not rotational invariant. For the real pro- 
cesses, SM"^ = 0, the second term vanishes for both channels, and the effective interaction 
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is independent on the cutoff function and coincide with the result of Tamm-Dancoff ap- 
proximation 

^eff - "i^^' (28) 

and similarly in annihilation channel. The same holds to the leading order of nonrela- 
tivistic expansion •p'^/m? <C 1 [15 . 



Estimate current-current term in exchange channel, Eq. which defines nominator 
of the interaction Eq. (]28|). We work in the Lepage-Brodsky convention for the spinors 

2 

A) = — (p^ + /?m a_LRL)A+XA 

,{E + j3m + ap)A+xx, (29) 



N 

and similarly for f (p, A) with the change m —>■ — m and x\ ~^ X~\ the above formula. 
The second expression for spinor holds quite in general for the solution of Dirac equation, 
where = E = y/p'^^frnF and the normalization factor = E + p^. Here /3 = 7", 
a = 7°7, A+ = 1/2(1 -|- a^) is the projection operator [|l|, and the spinor 



Xx 



. , 



is defined through the usual two-component spinors 

^^=(0) ' ^^(1)' 

Using the explicit representation for 7 matrices and projection operators A+,A_ and 
relations between them |T[, one has 

uip,Xh%ip',X') = -^^tiiE + p,)iE' + p[) + ipp')+t[pxp']a 

+ i[p X a]^{p'^ + m) -i[p' X a]^{p^ + m) + - Pzp'^)Cx' 
u{p, XhW, A') = -^^+{{E + pz){p'^ + i[p' X af) + {E' + p'J(p^ - i[p x a]') 

+ 5'\EE' -m^ -{pp')-i[pxp']a 

+ i[p X a]\E' - m) - i[p' X a]\E - m)) 

+ z£^V((E + p,)(m + p;)-(E' + pg(m + p,)))eA', (30) 

where i = 1,2,3; p = {pzi^i) is the three momentum and = e*-'^, e^'^ = —e"^^ = 1. 
Introducing three-momentum transfer and its mean 

q = p' — p 

k = \{P' + P), (31) 
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where 4(g^) = 6M'^ = 2{E^ - E'^) with 6M'^ defined above, Eq. (H) is written 

u{p,Xh\{p',X') = -±=^+{{E + h){E' + h)+k'-qy4-t[qxk]a 

- i[qx a],{k, + m) + q,{{E - E')/2 + i[k x a],) + w? - k^^y 
u{p,XhW,X') = -^^+{{{E + E')/2 + k,){2k^ + z[qxar) 

+ {(E-E')-qM/'2 + ^[kxaf) 

+ 5'\EE' - - {P - f/A) + i[q x k]a 

- t[qx afHE + E')/2 - m) - i[k x afiE - E')) 

+ ie''a\{k, + m){E-E')+q,{{E + E')/2-mmy, (32) 

where no approximations are done so far. Excluding the overall normalization factor the 



first lines in Eq. (|32D for scalar and vector current terms contain rotationally invariant 
parts (except terms proportional to kz), which coincide with the corresponding expressions 
when making use of Bjorken-Drell convention for spinors [jl2|. 



Merkel et.al. |T^ showed, that as far as the energy is conserved, this part gives rise 
to familiar spin dependent forces. The rest terms in Eq. (^) are obviously not rota- 
tionally invariant, particularly when the spacial rotations are performed perpendicular 
to the z-axis. Expanding expression Eq. (|32D to the second order in |p|/m -C 1 and 
performing the unitary transformation in spin space, Brisudova et.al. obtained Breit- 
Fermi spin-spin and tensor interactions. It seems to be impossible to reproduce full set 
of Breit-Fermi terms from the second order effective interaction in the light-front gauge. 
Also it is complicated to cover rotational symmetry on the level of light-front effective 
Hamiltonians without additional approximations are done. We use directly the effective 
electron-positron interaction Eq. ( P5| ) for numerical calculations of positronium spectrum. 
We aim to get fine structure and to investigate rotational symmetry on the level of spec- 
trum. The impact of different cutoff functions is also considered. The results of these 
calculations are presented in the next section. 



Spinors as used by Bjorken-Drell are 



A) = -^{E + 13m + ap)x\ , (33) 
V N 



2 



where N = E + m and xx is defined in the main text. The correspoding expressions for current terms 
are 

u(p,A)7%(p',A') = -^^+{{E + m){E' +m) + P-qy4-i[qxk]a)^y 



1 

^/NN' 

+ {E~E'){q'/2 + t[kxaY))^y, (34) 



For the energy conserving process this expression was obtained in p 
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n 


Term 


Betpt 


Be 


Bl 


Bg 


Bs 


1 


I'^o 


1.118125 


1.049550 


1.101027 


1.026170 


0.920921 


2 




0.998125 


1.001010 


1.049700 


0.981969 


0.885347 


3 




0.268633 


0.260237 


0.266490 


0.260642 


0.242607 


4 




0.253633 


0.253804 


0.259506 


0.254765 


0.234312 


5 


2^Pi 


0.253633 


0.257969 


0.263056 


0.257664 


0.237611 


6 


2'Po 


0.261133 


0.267070 


0.273826 


0.266563 


0.243075 


7 


2^Pi 


0.255508 


0.259667 


0.265412 


0.260127 


0.238135 


8 


23P2 


0.251008 


0.255258 


0.260345 


0.255498 


0.236383 



Table 1: Binding coefficients, -B„ = 4(2 — M„)/q;^ {a = 0.3), for the lowest modes of the 
positronium spectrum at = for the equal time perturbation theory up to order 
{Betpt) compared to our calculations with exponential {Be), Gaussian {Bq) and sharp 
{Bs) cutoffs. Bq is obtained using only g^^, part of interaction; for Bq 'rji^rj^' term is 
included. Exchange channel is considered. 



n 


Term 


Be 




Bg 


Bs 


1 




1.049550 


1.101270 


1.026170 


0.920921 


2 




0.936800 


0.978018 


0.921847 


0.834004 


3 




0.260237 


0.266490 


0.260642 


0.242624 


4 




0.255292 


0.260383 


0.255615 


0.234338 


5 




0.257969 


0.263056 


0.257664 


0.236383 


6 


2'Po 


0.267090 


0.273847 


0.266626 


0.243075 


7 


2^Pi 


0.259667 


0.265412 


0.260127 


0.237611 


8 


2^P2 


0.245615 


0.250821 


0.247091 


0.230901 



Table 2: Binding coefficients, Bn = 4(2 — Mn)/a^ {a = 0.3), for the lowest modes of the 
positronium spectrum at = for our calculations with exponential {Be), Gaussian 
{Bg) and sharp {Bs) cutoffs. Bq includes 'rj^^r]^^ term in exchange channel; Bq does not. 
Exchange and annihilation channels are considered. 



n 


Term 


SBe 


5Bg 




SBs 




2 


13^1 


6.30 10"^ 


1.76 10" 


-3 


1.18 10 


-3 


4 


23^1 


8.40 10"^ 


1.77 10" 


-4 


9.0 10" 


5 


5 


2^Pi 


-1.30 10"^ 


-7.47 10 


-4 


-9.1 10" 


-5 


7 


2^Pi 


-4.08 10-^ 


-4.08 10 


-4 


1.4 10" 


-4 


8 


2^P2 


5 10-6 


-7.7 10- 


5 


4.15 10 


-4 



Table 3: Difference in the corresponding energy levels between = and Jz = 1 states 
for exponential {5 Be), Gaussian {6Bg) and sharp {SBs) cutoffs. Exchange is channel is 
considered. 
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n 


Term 




X D 








2 


l^^i 


-1.411 10-3 


-7.86 10" 


-4 


-1.65 10" 


-3 


4 


23^1 


-4.1 10-5 


-4.0 10- 


5 


-1.15 10' 


~4 


5 


2iPi 


-6.4 10-5 


-6.52 10- 


-4 


-4.60 10' 


~4 


7 


23Pi 


-4.69 10-^ 


-4.74 10" 


-4 


-1.40 10" 


-4 


8 


23P2 


-1.96 10"^ 


-1.36 10" 


-4 


-2.44 10" 


-4 



Table 4: Difference in the corresponding energy levels between Jz = and Jz = 1 states for 
exponential {5 Be)-, Gaussian {6Bg) and sharp {6Bs) cutoffs. Exchange and annihilation 
channels are considered. 

4 Mass spectrum of positronium 

We solve the integral equation (P, with interaction kernel given in Eq.(|^), for positronium 
spectrum numerically. Effective interaction with different choice of cutoffs is sumarized 
in Appendix A. 

In polar coordinates the light-front variables are {k±]x) = (/t_L, x); therefore the 
matrix elements of the effective interaction Eq.(^ depend on the angles and if', i.e. 
(x, K±, {p; Ai, A2|Kff|2^'5 i^'_L, A'^, A2). In order to introduce the spectroscopic notation for 
positronium mass spectrum we integrate out the angular degree of freedom, introducing 
a discrete quantum number = n, n & 7a (actually for the annihilation channel only 
|<7z| < 1 is possible), 

(x, k_l; J2, Ai, A2|KfT|x', J^, A'^, A2) 




where L.^ = Jz — S^; = ^^ + ^- and the states can be classified (strictly speaking only for 
rotationally invariant systems) according to their quantum numbers of total angular mo- 
mentum J, orbit angular momentum L, and total spin S. Definition of angular momentum 
operators in light-front dynamics is problematic because they include interactions. 

The matrix elements of the effective interaction before integrating over the angles, 
(x, K^, Ai, A2|Kfr|a;', v^'; A'^, A2), and after the integration inroducing the total mo- 
mentum, J2, (x, Jz, Ai, A2|Kfr|a;', J'z, X'l, A2) for different cutoff functions are given 
in the exchange and annihilation channels in Appendices B and C, respectively. 

Now we proceed to solve for the positronium spectrum in all sectors of Jz- For this pur- 
pose we formulate the light-front integral equation Eq. (^ in the form where the integral 
kernel is given by the effective interaction for the total momentum Jz, Eq. (|35D . After the 
change of variables Eq. (0) we parametrize p = {n±,Pz) = (a* sin 9 cos ip, fi sin 9 sin ip, cos 9) 
The Jacobian of this transformation Eq. (|19D is given 

Im^ + jj"^ sin^ 9 

^ 2 (m2 + /i2)3/2 



(36) 
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One obtaines then the integral equation 

{Ml - A{m^ + /i^))V^„(/x, cos^; J„ Ai, A2 

+ E [ df^' r 



+1 , ^,/i'2m2 + ^'2(l-cos2^') 
rfcost^ 



2 (m2 + ^'2)3/2 

X (/i,cos^; J„Ai,A2|Kfr|/i',cos^'; J^,A;,A^)^„(/i',cos^'; J^,A;,A'2) = 0. (37) 



The integration domain D, defined in Eq. (13), is given now hj fie [0; f ]. Neither nor 
5*2 are good quantum numbers; therefore we set = Jz — S^. 

The integral equation Eq. ( pT] ) is used to calculate positronium mass spectrum numer- 
ically. Note, that if one succeeds to integrate out the angular degrees of freedom for the 
effective interaction Eq. (|35D analytically, one has 2-dimensional integration in Eq. (|3^) 
instead of 3-dimensional one in the original integral equation (Hf) to perform numerically. 



We use the numerical code |TT[, worked out by Uwe Trittmann for the similar problem 
This code includes for the numerical integration the Gauss-Legendre algorithm 
(Gaussian quadratures). To improve the numerical convergence the technique of Coulomb 
counterterms is included. The problem has been solved for all components of the total 
angular momentum, J^. 

Positronium spectrum is mainly defined by the Coulomb singularity 

(38) 

which is an integrable one analytically and also, by use of technique of Coulomb countert- 
erms, numerically. In this region 6Q^ —>■ and the energy denominator —>■ q"^ Eq. (|27|), 
giving rise to the leading order Coulomb behavior for the effective interaction Eq. (|28|) , 
independent on the cutoff function. We use therefore standard Coulomb counterterms, 
introduced for the Coulomb problem Eq. (|28D [|10| , pl|] , in the case of all cutoffs. Basing 
also on the argument Eq. (pHj), we expect the same pattern of levels for different cutoffs, 
that we prove numerically to be true. 

Another important limiting case to study effective interaction Eq. (^), namely its 
exchange part, is the coUinear limit 

(39) 



that is special for light-front calculations. Because of Eq. ( |26D the variable ^ ~ q \ 
resulting for the 'rj^rj'^ part of effective interaction to be 

^e'(l-^'(0)), (40) 

which is finite in this limit. This is true for the regular cutoff functions, as in the case 
of exponential and gaussian cutoffs, where the derivative d^{0)/d^ is well defined. For 
sharp cutoff this condition is not fulfilled, and the effective interaction contains the 1/q'^ 
type of singularity in this case (see Appendix A). We do not associate any physics with 
this singularity, considering it as a consequence of artificial choice of cutoff, which corre- 
sponds to singular generator of unitary transformation Eq. (|^). We omit the 'rj^rjl term 
in exchange channel for sharp cutoff in numerical calculations. 
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We argued that the region of Coulomb singularity, and hence 'g'^^^ part of effective 
interaction, determines mainly the positronium spectrum. However, including 'rjurjl part 
for gaussian cutoff shifts all levels as a whole of about 5 — 7%, since this part is diagonal in 
spin space (Appendix B), and improves the data to be near the result obtained in covariant 
equal time calculations (Table 1). Presumably, it is necessary to take into account 'r/^r/^ 
term in exchange channel also for sharp cutoff after the proper regularization of infrared 
longitudinal divergences is done. 

We place the results of calculations for three different cutoffs, performed in exchange 
and including both exchange and annihilation channels, in Tables 1 and 2, respectively. 
The corresponding set of figures is presented in Fig.l and Fig.2. We get the ionization 
threshold at ~ 4m^, the Bohr spectrum, and the fine structure. Including annihilation 
part increases the splittings twice as large for the lowest multiplets. 

As one can see from presents figures, certain mass eigenvalues at = arc degenerate 
with certain eigenvalues at other to a very high degree of numerical precision. As an 
example, consider the second lowest eigenvalue for = 0. It is degenerate with the lowest 
eigenvalue for = ±1, and can thus be classified as a member of the triplet with J — 1. 
Correspondingly, the lowest eigenvalue for Jz — having no companion can be classified 
as the singlet state with J = 0. Quite in general one can interpret degenerate multiplets 
as members of a state with total angular momentum J = 2Jz^max + 1- One can get the 
quantum number of total angular momentum J from the number of degenerate states for 
a fixed eigenvalue M^. One can make contact with the conventional classification scheme 
^^'^^Lj% as indicated in Tables 1-2. 

Such pattern of spectrum is driven by rotational invariance. To trace rotational sym- 
metry we calculate the difference of energy levels between Jz = and = 1 states for the 
lowest multiplets. The data are given for exchange and including annihilation channnel 
in Tables 3 and 4, respectively. Annihilation part makes corresponding states practically 
degenerate (see Tables 4 and Figure 2). 

5 Conclusion 

The numerical solution of positronium bound state problem, with the effective electron- 
positron interaction obtained by the flow equations, is presented. No approximations 
along numerical procedure are done. 

Concerning the spin-splittings the best agreement with covariant calculations is ob- 
tained for gaussian cutoff, the worst results arc for sharp cutoff. Rotational invariance 
is traced on the level of spectrum by studing the degree of degeneracy of corresponding 
states with the same total momentum but different projection Jz in the multiplet. Again, 
better results are obtained for exponential and gaussian cutoff functions than for sharp 
cutoff. This suggests, that smooth cutoff functions are preferable to perform calculations. 
Including annihilation channel improves the extend of degeneracy. 

For the sharp cutoff the lowest multiplet is placed higher than the one in case of 
exponential and gaussian cutoffs. The reason is in disregarding the infrared divergent 
part, which is diagonal in spin space and shifts the spectrum as a whole down. The 
question how to regularize this part and include it in mass spectrum calculations should 
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be considered. Generally, the impact of the different choice of cutoff functions on the 
spectrum is small. 

In this work wc solve the bound state integral equation for the one fixed integration in- 
terval. Integration domain introduces the ultraviolet- cutoff dependence of invariant mass 
squared M^(A), that reflects renormalization group properties of the effective coupling 
constant. We leave this question for the future study. 
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Figure 1: The invariant mass- squared spectrum for positronium versus tlie projection 
of the total spin, J^, excluding annihilation with exponential, Gaussian and sharp cutoffs. 
The number of integration points is A^i = A^2 = 21. 
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Figure 2: The invariant mass-squared spectrum for positronium versus the projection 
of the total spin, J^, including annihilation with exponential, Gaussian and sharp cutoffs. 
The number of integration points is A^i = A^2 = 21. 
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A Defining different cut-ofFs 

In this appendix we summarize the results for the effective electron-positron interaction, 
generated by the flow equations with different similarity functions. In the practical work, 
three different similarity function will be studied explicitly: 

(1) the exponential cut-off, (2) the gaussian cut-off, and (3) the sharp cut-off. 
(1) Exponential cut-off 



/(AO 
e{De,D,) 



exp 



\D\l) 



D. 



47]-2 



q+D 
1 



a 

47]-2 



p+D 
1 

M2 ' 



(41) 



where D = l/2{De + Dg) and D = l/2{Da + D^). The first choice of similarity function 
gives exactly the result of perturbation theory. 
(2) Gaussian cut-off 



/(AO 
e(A,A-) 



exp 



-D'l 



Dl 



D^^ + Dl 



2^ 



eff 



a 

4^ 

a 

4^ 

a 

4^ 

a 

47]-2 



1+e 

Qf^y De + De 



VfiVu {De - De 



Dl 



q+ Dl + Dl 2g+2 Dl 
p+Dl + Dl 2p+2 Dl + Dl 



Q\ 

M2 



p-1 



M4 



M4 + 



(42) 



where we understand under = {Q'^Y and 5Q^ = {6Q'^y with Q'^ and 6Q'^ defined in 
Eq. (P). 



(3) Sharp cut-off 

f{D;l) 
e{De,D,) 



^^(1- 
a 



\D\l) 

- \D,\) ; ^(0 = sign(0 



eff 



47]-2 



2q 



2 (De - D 



5W 



\De\-\De\) , e{\De\-\De\Y 



De 

\De\ - \De\) 



D-e 



\D, 



\De\r 



De 



D,. 



a 

47]-2 



IAD 



I A 



IA|)~ 



D, 



D. 
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{Da - D. 



a 
47r2 

a 



(7^70. 
(7"7''). 




\Df, 



\Da\y 



Df, 



Q2 + ^+2 n2 



M2 



M2 



Q2 + |5Q2| 

M2 + \SM^\ 



(43) 



The motivation to choose these cutoff functions is the following. Using exponential cutoff 
in flow equations one generates the same interaction as obtained also in Tamm-Dancoff 
approach, where numerical calculations of positronium spectrum are performed [1^, and 
we use this numerical code here. Note also, that for this cutoff the effective interaction 
looks very much as in covariant calculations: it contains only '(yf^,^' part, and 'rj^rj^' part 
is identically zero, so that there is no collinear problem. Gaussian cutoff corresponds to 
the original choice of generator Eq. (|^) by Wegner as commutator of diagonal, particle 
number conserving, and off-diagonal, particle number changing, parts of Hamiltonian. 
Sharp cutoff is used often in the alternative similarity scheme to perform calculations M . 
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B The matrix elements in the exchange channel 

In this Appendix we follow the scheme of the work |TD| to calculate the matrix elements 
of the effective interaction in the exchange channel.^ Here, we list the general, angle- 
dependent matrix elements defining the effective interaction in the exchange channel and 
the corresponding matrix elements of the effective interaction for arbitrary Jz, after in- 
tegrating out the angles. Exchange part of the effective interaction for three different 
cut-offs Eqs. (^l|-^3D can be written 



where explicitly one has 
(1) Exponential cut-off 



(2) Gaussian cut-off 



V,^ = -£^{YY)B,., (44) 



where Re and Im are real and imaginary parts, respectively, and = —1. 
(3) Sharp cut-off 



Q2-SQ2 Q2 + SQ2J 

where q = p[ —pi is the momentum transfer; and (7^7^^) for the exchange channel is given 
in Eq. @. We omit index 'ex' everywhere. 

It is convenient to extract the angular dependence in the functions 

QI = ai — boost 
Ql = (12 — b cos t 
t = (48) 

where we define 

k± = /c^ (cos 9?, sin 99) (49) 
in polar system; here the terms are given 



ai = — k, H k, + m 



X x' xx' 



^ Some of these calculations can be found in llif 
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= kl + k'l + ix-x')(kli-k-k'li-^)) ' -^^^ 

\ X X J 

a2 = ^ A;l + ^ -kl + ' 



l-x ^ l-x' ^ (1 -a;)(l -a;') 



^^^^ X -^l-xV ' ''^ (l-x)(l-x') 
h = 2k±k'^ (50) 

Then the functions in Eqs. (^^- |47|) are given 

= a ~ b cos t 

6Q^ = 6a, (51) 

where 

a = ^(ai + a2) 

5a = ^(01-02), (52) 

It is useful to display the matrix elements of the effective interaction in the form of 
tables. The matrix elements depend on the one hand on the momenta of the electron 
and positron, respectively, and on the other hand on their helicities before and after 
the interaction. The dependence on the helicities occur during the calculation of these 
functions E{x, Ai, A2IX', k'j_; X[, X'2) in part I and G{x, k±_; Ai, A2IX', /c^; X'^, X'2) in part 
II as different Kronecker deltas These functions are displayed in the form of helicity 
tables. We use the following notation for the elements of the tables 

F,(l,2) Ei{xM\x'fy- Gi{x,k^-x\k'^) (53) 

Also we have used in both cases for the permutation of particle and anti-particle 

F;{x, k^] x\ k'^) = ^3(1 - X, -k±; 1 - x', -k'^) (54) 

one has the corresponding for the elements of arbitrary J^; in the case when the function 
additionally depends on the component of the total angular momentum = n we have 
introduced 

Fi{n) = Fi{-n) (55) 



B.l The helicity table 

To calculate the matrix elements of the effective interaction in the exchange channel we 
use the matrix elements of the Dirac spinors listed in Table 1 ||^. Also the following holds 
V\'{p)Tvx{q) =ux{q)l°'ux'{p). 

We introduce for the matrix elements entering in the effective interaction Eqs. (^5|-^) 

2EW(x,fcx;Ai,A2|x',^l;A;,A;) = (7 V)^?/.. , (56) 
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M 



r 



r 



7i 



7i 



Vk+k 



=u{k',X')Mu{k,X) 



24 



k+k'- 



k'^ 



1 1 

' k+ 



Table 5: Matrix elements of the Dirac spinors. 

with {rr)9,u = 1(7+7-) + ^(7-7+) - (7?) - (7I) and 

2E^^\x,kr,Xi,X2\x',k'^;X[,X'^) = {^YKVu-^ 

with {'~i^'~i^)rj^rjy = {•y'^'j^); where 

k±; Xi)'j^'u{x', k'j_; X[)) {v{l - x', -k'^; X'^) Y^i^ - A2)) 



7+2 



'xx' 



^i\-x)(\-x') 



These functions are displayed in Table 



final : initial 


(A'i,A'2) =TT 


(A'i,A'2) =Ti 


(A'i,A'2) =iT 


(A'i,A'2) =ii 


(Ai,A2) =TT 


Ei(l,2) 


E3*(l,2) 


^3(1,2) 





(Ai,A2) =n 


^3(2,1) 


i?2(l,2) 


^4(1,2) 


-i?3(2,l) 


(Ai,A2) =iT 


i^3(2,l) 


^4(1,2) 


i^2(l,2) 


-^3(2,1) 


(Ai,A2) =ii 





-^3(1,2) 


-E3*(l,2) 


i^i(l,2) 



(57) 



(5^ 



Table 6: General helicity table defining the effective interaction in the exchange channel. 



The matrix elements i?|"''(l, 2) = E'f^\x^ k±; x', k'^j with n = 1 and = 2 for '(7' and 



rj^rji, terms, respectively, are the following 



E^^\x,kr,x',k'^) 



1 



k I k\ 



m 



XX' 



(1 — a;)(l — x') J xx'{l — x)(l — x') 
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E2 ^ (x, k±; X , 



m 



^xx' (1 — a;)(l — x') 



k'l- 



— x) x'{l — x') 



+ k±k\ I — + 



XX 



(1 -x') 



Ei'\x,kr,x\k'^) 
E^ ''(x, k±] X , 



m 
2 



1 - x' . 



, . k^e'^ - k^ ^ 
' ' 1 — X 

/\2 



(x — x') 



xx'(l — x)(l — x') 



(59) 



and 



El '{x,k±;x' 
Ef^(x,A?x;x',A?l) = Ei^\x,kr,x\k'^) =0 

B.2 The helicity table for arbitrary J^. 



(x — x')^ 



(60) 



Following the description given in the main text Eq. (^) we integrate out the angles 
in the effective interaction in the exchange channel. For the matrix elements of the 
effective interaction for an arbitrary = n with n G Z we introduce the functions 
G(x, k±; Ai, A2IX', k'j_; X[, Ag) = (x, k^; J^, Xi, A2|Kfr|a;', k'^; J'^, X[, A2) in the exchange chan- 
nel and obtain the helicity Table ^ 



final : initial 


(A'i,A'2) =TT 


(Ai,A',) =n 


(A'i,A'2) =iT 


(A;,A'2) =ii 


(Ai,A2) =TT 


Gi(l,2) 


G*(l,2) 


G3(l,2) 





(Ai,A2) =Ti 


G^(2,l) 


G2(l,2) 


^4(1, 2) 


-^3(2,1) 


(Ai,A2) =iT 


^3(2,1) 


G4(l,2) 


^2(1,2) 


-G^(2,l) 


(Ai,A2) =ii 





-^3(1, 2) 


-G*(l,2) 


Gi(l,2) 



Table 7: Helicity table of the effective interaction for Jz = ±n, x > x'. 
Here, the functions Gi{l, 2) = Gi{x, k±; x', k'j_) are given 

Gi{x,k±;x',k'j_) = 



^+ (l-x)(l-x0 ) 



+ 



ki k , 



26a 



xx'(l — x)(l 



— x'j (x — X 



G'2(x, k±; x', k'j_ 



k'l 



') J x(l — x) x'[l — x')J 



\ \xx' (1 — x)(l — x') 
+ *x*l - n|) + + «l)) 
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—Int{\n\) 



m ( , I , 1 — x' 



Gslx, k±; x' , kj) = I k^Int{\l — n\) — k± Int{\n\) 

xx' \ 1 — X ^ 

G^{x,k±;x',k'^_) = -^'^^7(YZ-^^r^^^^^(l^l) 



we define 



a f'^'^ cos nt 
I[n;a,b) = ——-^ / dt- 



2tt^ Jo a — b cos t 

then in Eq. (^) the following functions are introduced 
(1) Exponential cut-off 

Int{n) = I{n;a,b) 
Int{n) = 0, 



(2) Gaussian cut-off 



Int{n) = Iiel{n; a + i6a, b) 
Int{n) = lml{n; a + iSa, b) 



(3) Sharp cut-off 



Int{n) = 6{—6a)I{n; a — 6a,b) + 6{6a)I{n; a + 6a,b) 
Int{n) = 6{—6a)I{n; a — Sa,b) — 6{6a)I{n; a + Sa,b) 



also a + 5a = ai anda — 6a = 02- 
Explicitly is used 



, cos 1 fa — — 6^ \ 

dt = 2ti ^ 

a — boost — b'^ \ b J 

sin nt 

dt = , 

a — b cos t 

where a can contain imaginary part as in the case of gaussian cutoff. 
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C The matrix elements in the annihilation channel 

We repeat the same calculations for the matrix elements of the effective interaction in the 
annihilation channel. Annihilation part of the effective interaction can be written 

Ve^ = -£^{7^Y)C,., (67) 



where one has 



^ _ ± (<S)ab <S)ba\ V^^Vv 



a 



M2 1-/52 p+2 ' ^^^^ 

in the frame p± = 0.0 Explicitly the annihilation part of the effective interaction for 
different cut-offs Eq. (PI)- Eq. (^) is given 
(1) Exponential cut-off 



(2) Gaussian cut-off 



(3) Sharp cut-off 



where p'^ = pf + P2 is the total momentum; and {'j^'j'^) for annihilation is defined in 
Eq. (H). The functions present in Eq. (^9])- Eq. (^1]) are given in the light-front frame 



x'{l — x') 



2 I + 



= 7, (72) 

X 1 — X 



Indeed {l^'Y)9tJ^y = i(7^7") + + {l''Y)9i^u'^ therefore it holds 

_ ± r^t^jPy - Pu) + ^l^{P^i - Pjl) _ P^ 

gfif — 9au + j_ ,2 '?A''7i' J 

The 4- momentum of the photon in the t-channel can be written = + P2p ~ = pi^ + 

P2fj. — Vfj.Db/'^' with Da,Di, defined in Eq. dla). The Dirac equation {pi + p2)^u{pi)^^v{p2) — allows 
then to write Pfj,u{pi, Xi)j'^v{p2, X2) — — A'/^(2p+)7]^u(pi, Ai)7^w(p2, A2). Thus, when = 0, one has 

9^i^ 9^u -2^ ' 

p+ 

where the arrow means that this tensor should be contracted with (7^7'') in the annihilation channel. 
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we remind also 

6M' = \{Ml-Ml) (73) 

Note that the energy denominators of the effective interaction in the annihilation 
channel do not depend on the angles if, if' . 



M 


-jl=v{k\\')Mu{kA) 


7+ 




7~ 




5\, - m\ ( 




7i 


[k'+ ^ k+ ) 


5\, - mX ( 


,k'+ ^ k+) 




7i 




{l7+^k+, 


)4 



Table 8: Matrix elements of the Dirac spinors. 



C.l The helicity table 

For the calculation of matrix elements of effective interaction in the annihilation channel 
we use the matrix elements of the Dirac spinors listed in Table 4 0. Also the following 
holds {vx'iph°'ux{q)y = Ux{q)'J°'vx'{p)- 
We introduce 

2i/(i)(x, k^; Ai, X2\x', fcl; A^, A'^) = {YY)9i. = '{iD ' ill) 
2H^'\x,kr,XuX2\x',k'^;X[,X',) = {^YKv.^ (74) 

where 



^^i^Y) = (^(1 ~ >^'2h^u{x', k'^; X[)) {u{x, k±; Ai)7'"t;(l - x, -k±; A2) ^^^^ 

\Jx'{l — x') ^a;(l — x) 
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finalrinitial 


(a;,a'2) =tt 


iK,K) =n 


(a;,a'2) =iT 


(A'i,A^) =ii 


(Ai,A2) =TT 


Hi{l,2) 




h;{2,i) 





(Ai,A2) =n 


^3(1,2) 


H*{1,2) 


7/4(2, 1) 





(Ai,A2) =iT 




H,{1,2) 


i^2(l,2) 





(Ai,A2) =U 















Table 9: General helicity table defining the effective interaction in the annihilation chan- 
nel. 



These functions are displayed in the Table 5. 

Here, the matrix elements h\^\i^ 2) = hI^\x, k±; x', k'j_) are the following 

Hi'\x,kr,x'X) = -^'(- + 7^) (A + 7^, 

\x 1 — xj \x' 1 — x' 

(1) ^ / ^' / e^if-f' 



\x 1 — xJ 1 — x' 
Hi'\x,kr,x',k^) = k^k'^ I ^^^^ _ J (76) 



and 



Hf\x, kr, X, k'^) = Hf\x, k±; x', fc^) = 

H^^\x, k±; x', k'^) = Hi^\x, kr, x', k'^) = 2 (77) 

C.2 The helicity table for | < 1 

The matrix elements of the effective interaction for Jz > F(x, A;^; Ai, A2|x', A;^; A'^, A2) = 
(,T, fcj^; J2, Ai, \2\Veff\x' ^ k'^; J'^, \[, A2) in the annihilation channel (the sum of the gener- 
ated interaction for = +1 and instantaneous graph for = 0) are given in Table 
6. 

The function Fj(l, 2) — Fi{x,k±; x' , k'j) are the following 

F,{x,kr,x,k,) = 

F,ix,kr,x',k',) = ^(^^^k.l.i + 25..0) 

F,ix,kr,x,k^) = 
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finalrinitial 


=TT 


(a;,a^) =Ti 


(a;,a'2) =iT 


(A'i,A^) =ii 


(Ai,A2) =TT 


^i(l,2) 


i^3(2,l) 


F3*(2,l) 





(Ai,A2) =n 


i^3(l,2) 


F*{1,2) 


i^4(2,l) 





(Ai,A2) =iT 


F3*(l,2) 


i^4(l,2) 


i^2(l,2) 





(Ai,A2) =U 















Table 10: Helicity table of the effective interaction in the annihilation channel for > 0. 



/ a f 1 k±k' ^ 



where we have introduced 
(1) Exponential cut-ofF 



(2) Gaussian cut-off 



11 



1 _ M'^+M^ 
Q, Mf; + Mt ' 



(3) Shcirp cut-ofF 



Q Ml M, 



(78) 



(79) 



(80) 



1 _ e{Mi - M^) e{M^ - Ml) _ 



The table for = —1 is obtained by inverting all helicities, i.e. 

F{J, = +1; Ai, A2) = -AiF(J, = -1; -Ai, -A2) , (82) 

The matrix elements of the effective interaction in the annihilation channel are nonzero 
only for | J^l < 1 due to the restriction on the angular momentum of the photon. 
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